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We study the effect of turbulent drift of a large-scale magnetic field that results from the interaction of helical convective motions and 
differential rotation in the solar convection zone. The principal direction of the drift corresponds to the direction of the large-scale vorticity 
vector. Thus, the effect produces a latitudinal transport of the large-scale magnetic field in the convective zone wherever the angular 
velocity has a strong radial gradient. The direction of the drift depends on the sign of helicity and it is defined by the Parker-Yoshimura 
rule. The analytic calculations are done within the framework of mean-field magnetohydrodynamics using the minimal r-approximation. 
We estimate the magnitude of the drift velocity and find that it can be several m/s near the base of the solar convection zone. The 
implications of this effect for the solar dynamo are illustrated on the basis of an axisymmetric mean-field dynamo model with a subsurface 
shear layer. We find that the helicity-vorticity pumping effect can have an influence on the features of the sunspot time-latitude diagram, 
producing a fast drift of the sunspot activity maximum at the rise phase of the cycle and a slow drift at the decay phase of the cycle. 
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1 Introduction 

It is believed that the evolution of the large-scale magnetic field of the Sun is governed by the interplay 
between large-scale motions, like differential rotation and meridional circulation, turbulent convection flows 
and magnetic fields. One of the most important issues in solar dynamo theory is related to the origin of 
the equatorial drift of sunspot activity in the equatorial regions and, simultaneously at high latitudes, 



the poleward drift of the location of large-scale unipolar regions and quiet prominences. Parker (1955) 



and Yoshimura (1975) suggested that the evolution of large-scale magnetic activity of the Sun can be 



interpreted as dynamo waves propagating perpendicular to the direction of shear from the differential 
rotation. They found that the propagation can be considered as a diffusion process, which follows the 
iso-rotation surfaces of angular velocity in the solar convection zone. The direction of propagation can be 
modified by meridional circulation, anisotropic diffusion and the effects of turbulent pumping (see, e.g., 
Choudhuri et al.p.995 , Kitchatinovj|2002 , Guerrero and de Gouveia Dal Pino||2008 ). The latter induces an 
effective drift of the large-scale magnetic field even though the mean flow of the turbulent medium may 
be zero. 

The turbulent pumping effects can be equally important both for dynamos without meridional circulation 
and for the meridional circulation-dominated dynamo regimes. For the latter case the velocity of turbulent 
pumping has to be comparable to the meridional circulation speed. It is known that an effect of this 
magnitude can be produced by diamagnetic pumping and perhaps by so-called topological pumping. Both 
effects produce pumping in the radial direction and have not a direct impact on the latitudinal drift of the 
large-scale magnetic field. 



Recently rtPipin 2008, Mitra et al. 2009, Leprovost and Kim 2010), it has been found that the helical 



convective motions and the helical turbulent magnetic fields interacting with large-scale magnetic fields 
and differential rotation can produce effective pumping in the direction of the large-scale vorticity vector. 
Thus, the effect produces a latitudinal transport of the large-scale magnetic field in the convective zone 
wherever the angular velocity has a strong radial gradient. It is believed that these regions, namely the 
tachocline beneath the solar convection zone and the subsurface shear layer, are important for the solar 
dynamo. Figure [1] illustrates the principal processes that induce the helicity-vorticity pumping effect. It is 
suggested that this effect produces an anisotropic drift of the large-scale magnetic field, which means that 
the different components of the large-scale magnetic field drift in different directions. Earlier work, e.g. 



by Kichatinov (1991) and Kleeorin and Rogachevskii (2003), suggests that the effect of anisotropy in the 
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Figure 1. The field lines of the large-scale magnetic field, B^ T \ are transformed by the helical motions to a twisted i?-like shape. 
This loop is folded by the large-scale shear, V^ T \ into the direction of the background large-scale mag netic field, B^ T \ The induced 
electromotive force has a component, S^ p \ which is perpendicular to the field B^ T \ The resulting effect is identical to the effective 
drift of the large-scale magnetic field along the x-axis, in the direction opposite to the large-scale vorticity vector W = V x V^ T \ i.e., 
£(P) ~ -W x B( T ). 



transport of mean-field is related to nonlinear effects of the global Coriolis force on the convection. Also, 



nonlinear effects of the large-scale magnetic field result in an anisotropy of turbulent pumping (Kleeorin 



et al. 1996). It is noteworthy, that the helicity- vorticity effect produces an anisotropy of the large-scale 



magnetic field drift already in the case of slow rotation and a weak magnetic field. A comprehensive study 
of the linear helicity- vorticity pumping effect for the case of weak shear and slow rotation was given by 



Rogachevskii et al. 



Brandenburg et al. 



2011) and their results were extended by DNS with a more general test-field method 



(2012 



In this paper we analytically estimate the helicity-vorticity pumping effect taking into account the 
Coriolis force due to global rotation. The calculations were done within the framework of mean-field 
magnetohydrodynamics using the minimal r-approximation. The results are applied to mean field dynamo 
models, which are used to examine this effect on the dynamo. The paper is structured as follows. In the 
next section we briefly outline the basic equations and assumptions, and consider the results of calculations. 
Next, we apply the results to the solar dynamo. In Section 3 we summarize the main results of the paper. 
The details of analytical calculations are given in the Appendices A and B. 



2 Basic equations 



In the spirit of mean- field magnetohydrodynamics, we split the physical quantities of the turbulent con- 
ducting fluid into mean and fluctuating parts where the mean part is defined as an ensemble average. One 
assumes the validity of the Reynolds rules. The magnetic field B and the velocity V are decomposed as 
B — B + b and V — V + u, respectively. Hereafter, we use small letters for the fluctuating parts and 
capital letters with an overbar for mean fields. Angle brackets are used for ensemble averages of products. 
We use the two-scale approximation (jRoberts and Soward 1975, Krause and Radler 1980) and assume 



that mean fields vary over much larger scales (both in time and in space) than fluctuating fields. The 
average effect of MHD-turbulence on the large-scale magnetic field (LSMF) evolution is described by the 
mean-electromotive force (MEMF), £ = (u X b). The governing equations for fluctuating magnetic field 
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and velocity are written in a rotating coordinate system as follows: 

db 

- = Vx(wxB + yx5)+ r]V 2 b + 0, (1) 

at v y 

d ^+2(f2xu) i = -V i ^p+^p-^j +vA Ui (2) 

+ - Vj (Bjbi + Bibj) - Vj (VjUi + ViUj) +fi + &, 
ji 

where stand for nonlinear contributions to the fluctuating fields, p is the fluctuating pressure, ft 
is the angular velocity responsible for the Coriolis force, V is mean flow which is a weakly variable in 
space, and / is the random force driving the turbulence. Equations and ^ are used to compute the 
mean-electromotive force £ = (u X b). It was computed with the help of the equations for the second 
moments of fluctuating velocity and magnetic fields using the double-scale Fourier transformation and the 
minimal r-approximations and for a given model of background turbulence. To simplify the estimation 
of nonlinear effects due to global rotation, we use scale-independent background turbulence spectra and 
correlation time. Details of the calculations are given in Appendix A. In what follows we discuss only those 
parts of the mean-electromotive force which are related to shear and the pumping effect. 



2.1 Results 

The large-scale shear flow is described by the tensor Vij = VVV It can be decomposed into a sum of strain 
and vorticity tensors, V jVi = - (Vij + Vj : i) — -£ij P W p , where W = V x V is the large-scale vorticity 

Z Z 

vector. The joint effect of large-scale shear, helical turbulent flows and magnetic fields can be expressed 
by the following contributions to the mean-electromotive force (omitting the a-effect): 

= (WxB) (fMhc + f^hjc) rl + V (B) (f^hc + f^h K ) rl (3) 
+ e [(e x W) B] (f^hc + f^h,c) rl + (e x W) (e-B) (f^h c + f^h K ) r c 2 , 

where V (B) = — - {Vij + Vj^) 1 e = — — is the unit vector along the rotation axis, r c is the typical 

relaxation time of turbulent flows and magnetic fields, Ifi^ — {u^^ X and = - 

are kinetic and current helicity of the background turbulence. These parameters are assumed to be known 
in advance. Functions fn(f2*) are given in Appendix B, they depend on the Coriolis number i?* = 2Q^r c 
and describe the nonlinear effect due the Coriolis force, and is the global rotation rate. 
For slow rotation, f2* <C 1, we perform a Taylor expansion of f&\f2*) and obtain 

£W = T l(WxB) (h c - h K ) + V (B) (3% - y hA . (4) 



The coefficients in the kinetic part of Eq. Q are two times larger than those found by |Rogachevskii 



et al. (2011). This difference results from our assumption that the background turbulence spectra and the 
correlation time are scale-independent. The results for the magnetic part are in agreement with our earlier 
findings (see |Pipin||2008 ). The first term in Eq. ^ describes turbulent pumping with an effective velocity 



r 2 W 

- £ — — (he — hjc) and the second term describes anisotropic turbulent pumping. Its structure depends on 
z 

the geometry of the shear flow. For large Coriolis numbers, J?* ^> 1, only the kinetic helicity contributions 
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Figure 2. The dependence of the pumping effects on the Coriolis number. Solid lines show contributions from kinetic helicity and dashed 
lines the same for current helicity. 



survive: 



2 2 
6 v 7 5 



(5) 



Figure |2]j$how the dependence of the pumping effects on the Coriolis number. We observe that for the 
terms (W x S) and V (B) the effects of kinetic helicity are non-monotonic and have a maximum at 
J?* « 1. The effects of current helicity for these terms are monotonically quenched with increasing values 
of J?*. The additional contributions in Eq. (3]) are rather small in comparison with the main terms. Thus, 
we can conclude that the first line in Eq. ([3) describes the leading effect of pumping due to the helicity of 
turbulent flows and magnetic field. Below, we drop the contributions from the second line in Eq. ^ from 
our analysis. 



2.2 Helicity -vorticity pumping in the solar convection zone 



2.2.1 The dynamo model To estimate the impact of this pumping effect on the dynamo we consider 
the example of a dynamo model which takes into account contributions of the mean electromotive force 
given by Eq. ([3]). T he dynamo model employed in this paper has been described in detail by Pipin and 
Kosovichev ( 2011a|c ). This type of dynamo was proposed originally by Brandenburg (2005). The reader 
may find the discussion for different types of mean- field dynamos in Brandenburg and Subramanian (2005) 
andlTobias and Weissl ([2007| . 



We study the standard mean-field induction equation in a perfectly conducting medium: 



dB 

~dt 



V x (S + U x B) 



(6) 



where S — u X b is the mean electromotive force, with u, b being fluctuating velocity and magnetic field, 
respectively, U is the mean velocity (differential rotation and meridional circulation), and the axisymmetric 
magnetic field is: 



B = e 6 B + V 



Ae d 



r sm( 



where 6 is the polar angle. The expression for the mean electromotive force S is given by|Pipinl(p008l). It 
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is expressed as follows: 



Si 



a 



13 



+ i\f> ) B - rjijkVjBk + S, 



(7) 



The new addition due to helicity and mean vorticity effects is marked by S . The tensor a^j represents 
the a-effect. It includes hydrodynamic and magnetic helicity contributions, 



<Mj = C a sin 6a\- + a 



2 nAH) . AM) 



13 



X - ■ 

13 



^ ' = % (f$ (e-AW) + /<?> (cAW)) } + 
+ e*i{3*r(/f° (e.^)) + /i°) (,AW))} + 



(8) 
(9) 
(10) 

(11) 



where the hydrodynamic part of the a-effect is defined by ol^\ A^ = Vlogp quantifies the density 
stratification, = Vlog [Vtj quantifies the turbulent diffusivity variation, and e = ft I \ft\ is a unit 



vector along the axis of rotation. The turbulent pumping, ")\f \ depends on mean density and turbulent 
diffusivity stratification, and on the Coriolis number f2* = 2r c i?o where r r is th e typical convective turnover 
time and f2o is the global angular velocity. Following the results of 



Pipin 



(2008), 7-^ is expressed as follows: 



A) 



3 VT {fi a) A^ + f[ a) (e-AW) e n ) e inj - ^ T f[ a) e 3 e inm e n At\ 



-in j- 



(12) 
(13) 



- 3 VT (e - 1) {i a) A^ + /« (e.iW)e n } 
The effect of turbulent diffusivity, which is anisotropic due to the Coriolis force, is given by: 

Vijk = 3t?t { {zf[ a) - s ijk - 2f[ a) eie n e njk + eC^f^ ej8 ik ^ . (14) 



14) describes Radler's ft x J effect. The functions f^-\\\ depend on the Coriolis 



(2011 



l i-ii) 



The last term in Eq. 

number. They can be found in Pipin (2008); see also Pipin and Kosovichev (2011a|) or Pipin and Sokoloff 



)). In the model, the parameter e = 



, which measures the ratio between magnetic and kinetic 



energies of the fluctuations in the background turbulence, is assumed to be equal to 1. This corresponds 



to perfect energy equipartition. The e contribution in the second line of Eq. (12) describes the param- 



agnetic effect (Kleeorin and Rogachevskii 2003). In the state of perfect energy equipartition the effect of 



diamagnetic pumping is compensated by the paramagnetic effect. We can, formally, skip the second line 



in Eq. (12) from our consideration if e — 1. To compare the magnitude of the helicity-vorticity pumping 



effect with the diamagnetic effect we will show results for the pumping velocity distribution with e — 0. 

The contribution of small-scale magnetic helicity x — a 'b (a is the fluctuating vector-potential of the 
magnetic field) to the o:-effect is defined as 



a 



(M) 



2/ 2 (a) * 



popi 2 



p p£ 2 ' 



(15) 



The nonlinear feedback of the large-scale magnetic field to the a-effect is described by a dynamical quench- 
ing due to the constraint of magnetic helicity conservation. The magnetic helicity, x > subject to a con- 
servation law, is described by the following equation (Kleeorin and Rogachevskii 1999, Subramanian and 
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Brandenburg||2004| ) : 



9 f t = -2 (S-B) - - JL. + V- (r? x Vx) , 



R x t c 



(16) 



where r c is a typical convective turnover time. The parameter R x controls the helicity dissipation rate 
without specifying the nature of the loss. The turnover time r c decreases from about 2 months at the 
bottom of the integration domain, which is located at O.71i? , to several hours at the top boundary 
located at O.99i? . It seems reasonable that the helicity dissipation is most efficient near the surface. The 



last term in Eq. (16) describes a turbulent diffusive flux of magnetic helicity (Mitra et al. 2010). 



We use the solar convection zone model computed by Stix (2002 ), in which the mixing-length is defined as 
£ = «mlt | , where = Vlogp quantifies the pressure variation, and q^mlt = 2. The turbulent 



diffusivity is parameterized in the form, t\t — C v r]^\ where f]^ J = is the characteristic mixing- 

o 

length turbulent diffusivity, i is the typical correlation length of the turbulence, and C n is a constant to 
control the efficiency of large-scale magnetic field dragging by the turbulent flow. Currently, this parameter 
cannot be introduced in the mean-field theory in a consistent way. In this paper we use — 0.05. The 
differential rotation profile, Q = ftofn m) (shown in Figj3^i) is a slightly modified version of the analytic 
approximation proposed by Antia et al.| ( 1998| ): 



,(°) 



fa (x, n) = -J- [f2 + 55 (x - 0.7) <f> (x, x ) </> (-x, -0.96) 
- 200 {x -0.95)0 (a;, 0.96)) 

+ (21P 3 (p) + 3P 5 (M)] (J^ + ^—tK) M, 



(17) 



jp (x) j e (x) 



l + exp(^)' Je l + exp(0fc)' 



where = 2.87 • 10 6 s 1 is the equatorial angular velocity of the Sun at the surface, x = r/i?©, 
(j) (x, x ) = 0.5 [1 + tanh [100 (x - x )]] : and x = 0.71. 



2.2.2 Pumping effects in the solar convection zone. The components of the strain tensor V in a 
spherical coordinate system are given by the matrix: 




where we take into account only the azimuthal component of the large-scale flow, Vjy,^) = r sin 6d r f2 (r, 0), 
V(Q,<p) = sinfldfli? (r, fl), so V (B) = (BV^ r ^ BV( ^ B^Vg Substituting this into Eq. ( 3) we find the 
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Figure 3. Distributions of the angular velocity and the turbulent parameters, and the kinetic helicity inside the solar convection zone. 
The bottom panel shows the patterns of t he pumpin g velocity fields for the toroidal magnetic fleld(left) and for the poloidal fleld(right). 
They were computed on the basis of Eqs. ( |18|19|2d] ). 

components of the mean-electromotive force for the helicity-vorticity pumping effect, 



s, 



£ 



(H) 



4 H) = ^ sin* <j [h K (/W - /«) + h c (/W - /W)] x ^-2 (i2 - l)[h K f™ + h c m B , (18) 



2 

Q*T, 



dx 



c J sin 2 ( 



i?*r r sin 9 



x 



hK (/i 7) - fl^+hc - / 2 (7) )]f^ - 2 M (tf - l)[%A W + Ac/W] b, (19) 



%(/ 3 (7) +/i (7) )+^ (/| 7) + ^ 7) 



d(n,A 

d (x, /i) 



(20) 



i? — 1 I j?*r r r 



xsint 



(7)' 



9 A sin 2 <9<9,4 N 



where /i£ — Cq 



X 



It remains to define the kinetic helicity distribution. We use a formula proposed in 



our earlier study (see |Kuzanyan et al.||2006 ), 



v (°) 2 d ( I \ 
hjc = C vCjc — ^; log ( pV^ (0)2 J *icos0, 



where F\ (J?*) was defined in the above cited paper. The radial profile of is shown in Figure 

cos# 



The 



radial profile of kinetic helicity is shown in Figure 3a of the above cited paper. The parameters Cy^c are 
introduced to switch on/off the pumping effects in the model. 
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Figure 4. The patterns of the total (including the diamagnetic and the density gradient effects) pumping velocity fields for the toroidal 
magnetic field(left) and for the poloidal field(right). 



The expressions given by Eq. O) are valid for the case of weak shear, when r c max (|ViVj|) C 1. In 



terms of the strain tensor V this condition of weak shear implies i?* max 



rd r Q 



Oaf? 



< 1. This is not 



valid at the bottom of the solar convect ion zone where the ra dial gradient of the angular velocity is strong 



and i?* ^> 1 and r c max (| ViVj |) « 2 

quenched with increasing shear inversely proportional to (r c max (| ViV^D) 
an ad- hoc quenching function for the pumping effect: 



Leprovost and Kim (2010) suggested that this pumping effect is 

1 2 

. Therefore, we introduce 



f(S) 



1 + C S Q* 



dQ 

dr 



+ 



dQ 



(21) 



where Cs is a constant to control the magnitude of the quenching, and s = 1. Results by |Leprovost and" 
Kim (2010) suggest 1 < s < 2 in relation to geometry of the large-scale shear. We find that for the solar 
convection zone the amplitude of the pumping effect does not change very much (~ 1 m/s) with s varying 
in the range 1 ... 2. 

From the given relations, using 

£ (H) = jjieS) x B 

we find the effective drift velocity, due to 

the helicity-vorticity pumping effect. Taking into account the variation of turbulence parameters in the 
solar convection zone we compute The bottom panel of Figure [3] shows the distribution of the 

velocity field JJ^ e ^^ for the helicity-vorticity pumping effect for the toroidal and poloidal components of 
the large-scale magnetic field. The maximum velocity drift occurs in the middle and at the bottom of the 
convection zone. The direction of drift has equatorial and polar cells corresponding to two regions in the 
solar convection zone with different signs of the radial gradient of the angular velocity. The anisotropy in 
transport of the toroidal and poloidal components of the large-scale magnetic field is clearly seen. 

The other important pumping effects are due to mean density and turbulence intensity gradients (Zel- 
dovich 1957, Kichatinov 1991, Kichatinov and Riidiger 1992, Tobias et al. 2001). These effects were esti- 
mated using Eq. (12). For these calculations we put C v = 1, e = 0, Cjc = 1 and % = 0. Figure |i] shows the 
sum of the pumping effects for the toroidal and poloidal components of mean magnetic fields including 
the helicity-vorticity pumping effect. In agreement with previous studies, it is found that the radial direc- 
tion is the principal direction of mean-field transport in the solar convection zone. In its upper part the 
transport is downward because of pumping due to the density gradient ( Kichatinov||1991 ). At the bottom 
of the convection zone the diamagnetic pumping effect produces downward transport as well (Kichatinov 



1991 , |Riidiger and Brandenburg 1995). The diamagnetic pumping is quenched inversely proportional to 
the Coriolis number (e.g., Kichatinov 1991, Pipin [2008]) and it has the same order of magnitude as the 



helicity-vorticity pumping effect. The latter effect modifies the direction of effective drift of the toroidal 
magnetic field near the bottom of the convection zone. There is also upward drift of the toroidal field at 
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dynamo period of the model. 
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is the maximum of the toroidal magnetic field strength inside the convection zone, P is the 
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low latitudes in the middle of the convection zone. It results from the combined effects of density gradient 



and global rotation ( jKichatinov 1991, Krivodubskij 2004). For the poloidal magnetic field the transport 
is downward everywhere in the convection zone. At the bottom of the convection zone the action of the 
diamagnetic pumping on the meridional component of the large-scale magnetic field is amplified due to 
the helicity-vorticity pumping effect. 

The obtained pattern of large-scale magnetic field drift in the solar convection zone does not take into 
account nonlinear effects, e.g., because of magnetic buoyancy. The effect of mean- field buoyancy is rather 
small compared with flux-tube buoyancy QKichatinov and Pipin|[l993 , cf. Guerrero and Kapyla||20lT ). 

To find out the current helicity counterpart of the pumping effect we analyze dynamo models by solving 

Eqs. jf}, 16). The governing parameters of the model are C n = 0.05, C u = \p a . We discuss the choice of 

the governing parameters later. The other parameters of the model are given in the Table [T] Because of 
the weakening factor C n the magnitude of the pumping velocity is about one order of magnitude smaller 
than what i s shown in Figure [4) 

Following Pipin and Kosovichev (2011b), we use a combination of "open" and "closed" boundary con- 
ditions at the top, controlled by a parameter 5 — 0.95, with 



5)£ e = 0. 



(22) 



This is similar to the boundary condition discussed by [Kitchatinov et al. (2000). For the poloidal field we 
apply a combination of the local condition A = and the requirement of a smooth transition from the 
internal poloidal field to the external potential (vacuum) field: 



Qj^ivac) 



dr 



+ (1-(J)A = 0. 



(23) 



We assume perfect conductivity at the bottom boundary with standard boundary conditions. For the 
magnetic helicity, similar to Guerrero et al. (2010), we put V r x = at the bottom of the domain and 
X = at the top of the convection zone. 
In this paper we study dynamo models which include Radler's 4? x J dynamo effect due to a large-scale 



current and global rotation (|Radler 



current (Rogachevskii and Kleeorin 



1969). There is also a dynamo effect due to large-scale shear and 
2003). The motivation to consider these addional turbulent sources 



in the mean- field dynamo comes from DNS dynamo experiments (Brandenburg and Kapyla 2007, Kapyla 



eFaLl[2008l|Hughes and ProctorpOO^ |Kapyla et al.||2009| ) and from our earlier studies (Pipin and Seehafer 
2009, Seehafer and Pipin||2009 ) . The dynamo effect due to large-scale current gives an additional source of 
large-scale poloidal magnetic field. This can help to solve the issue with the dynamo period being otherwise 
too short. Also, in the models the large-scale current dynamo effect produces less overlapping cycles than 
dynamo models with a-effect alone. The choice of parameters in the dynamo is justified by our previous 
studies ( Pipin and Seehafer|2009 , Pipin and Kosovichev|20lTc ), where we showed that solar-types dynamos 
can be obtained for C a /C u > 2. In those papers we find the approximate threshold to be C a « 0.02 for a 
given diffusivity dilution factor of C n — 0.05. 

As follows from the results given in FigjlJ the kinetic helicity-vorticity pumping effect has a negligible 
contribution in the near-surface layers, where downward pumping due to density stratification dominates. 
Therefore, it is expected that the surface dynamo waves are not affected if we discard magnetic helicity 
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Figure 5. The time- latitude diagrams for the toroidal and radial magnetic fields for the models Dl and D2: a) the model Dl, the toroidal 
field (iso-contours, ±.25KG) near the surface and the radial field (gray-scale density plot); b) the model Dl, the toroidal field at the 
bottom of the solar convection zone, the contours drawn in the range ±.5KG; c) the same as for item a) for the model D2; d) the same 
as for item b) for the model D2. 

from the dynamo equations. Figure [5] shows time-latitude diagrams for toroidal and radial magnetic fields 
at the surface and for toroidal magnetic field at the bottom of the convection zone for two dynamo models 
Dl and D2 with and without the helicity-vorticity pumping effect, but magnetic helicity is taken into 
account as the main dynamo quenching effect. To compare with observational data from a time-latitude 
diagram of sunspot area (e.g., Hathaway||2011 ), we multiply the toroidal field component B by factor sin#. 
This gives a quantity, which is proportional to the flux of large-scale toroidal field at colatitude 9. We 
further assume that the sunspot area is related to this flux. 

Near the surface, models Dl and D2 give similar patterns of magnetic field evolution. At the bottom 
of the convection zone model Dl shows both poleward and equatorward branches of the dynamo wave 
propagation that is in agreement with the Parker- Yoshimura rule. Both branches have nearly the same time 
scale that equals ^16 years. The results from model D2 show that at the bottom of the convection zone the 
poleward branch of the dynamo wave dominates. Thus we conclude that the helicity-vorticity pumping 
effect alters the propagation of the dynamo wave near the bottom of the solar convection zone. We find 
that models with magnetic helicity contributions to the pumping effect do not change this conclusion. 

Figure [6] shows a typical snapshot of the magnetic helicity distribution in the northern hemisphere for all 
our models. The helicity has a negative sign in the bulk of the solar convection zone. Regions with positive 
current helicity roughly correspond to domains of the negative large-scale current helicity concentration. 
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Figure 6. Snapshots for the mean magnetic field and the current helicity distributions at the north hemisphere in the model D4. Left 
panel shows the field lines of the poloidal component of the mean magnetic field. The right panel shows the toroidal magnetic field 
(iso-contours ±500G) and the current helicity (gray scale density plot). 



They are located in the middle of the solar convection zone and at the high and low latitudes near the 
top of the solar convection zone. As follows from Fig. |6j the pumping effect due to current helicity may 
be efficient in the upper part of the solar convection zone where it might intensify the equatorial drift of 
the dynamo wave along iso-surfaces of the angular velocity. We find that the pumping effect that results 
from magnetic helicity is rather small in our models. This may be due to the weakness of the magnetic 
field. Observations ( Zhang et al.||2010 ) give about one order magnitude larger current helicity than what 



X 

is shown in Fig. 6 In the model we estimate the current helicity as He = — 777. This result depends 

□ m Ho? 2 

essentially on the mixing length parameter £. The stronger helicity is concentrated to the surface, the 

larger He- In observations, we do not know from were the helical magnetic structures come from. In view 

of the given uncertainties we estimate the probable effect of a larger magnitude of magnetic helicity in 

the model by increasing the parameter Cq to 10 (model D3). In addition, we consider the results for the 

nonlinear model D4. It has a higher C a and a lower R x to increase the nonlinear impact of the magnetic 

helicity on the large-scale magnetic field evolution. 

The top panel of Figure [7] shows a time-latitude diagram of toroidal magnetic field and current helicity 
evolution near the surface for model D4. We find a positive sign of current helicity at the decay edges 
of the toroidal magnetic field butterfly diagram. There are also areas with positive magnetic helicity at 
high latitudes at the growing edges of the toroidal magnetic field butterfly diagram. The induced pumping 
velocity is about 1 cms -1 . The increase of the magnetic helicity pumping effect by a factor of 10 (model 
D3) shifts the latitude of the maximum of the toroidal magnetic field by about 5° toward the equator. The 
induced pumping velocity is about 5 cms -1 . 

Stronger nonlinearity (model D4) and a stronger magnetic helicity pumping effect (model D3) modify the 
butterfly diagram in different ways. Model D3 shows a simple shift of the maximum of toroidal magnetic 
field toward the equator. Model D4 shows a fast drift of large-scale toroidal field at the beginning of a 
cycle and a slow-down of the drift velocity as the cycle progresses. 

Figure [8] shows in more detail the latitudinal drift of the maximum of the toroidal magnetic field evolution 
during the cycle (left panel in the Figure [8]), 

A max (t) = 90° - max (\B s (9)\sm9) , (24) 

#>45° 



and the latitudinal drift of the centroid position of the toroidal magnetic field flux (cf. |Hathaway||2011 ) 



n fn 0Bs(0) sin 9d0 , . 

A C t = 90° - J \ /2 SK J , 25 

f* /2 B s {9) sin 0d0 

where B$(0) = (B (r, 0))( O 9 99) r * s ^ ne toroidal magnetic field, which is averaged over the surface layers. 
Note that the overlap between subsequent cycles influences the value of Ac more than the value of A max . 
The behaviour of A max in models D1,D2 and D3 reproduces qualitatively the exponential drift of maximum 



March 4, 2013 17:0 Geophysical and Astrophysical Fluid Dynamics pipin3 



12 



V. V. Pipin 




50 60 
TIME, [YR] 

Figure 7. Top, the near-surface time-latitude diagrams for the toroidal magnetic field and the current helicity for the models D4. Bottom, 
the near-surface time-latitude diagrams for the toroidal magnetic field and the latitudinal component of the drift velocity induced by the 
magnetic helicity for the model D3. 




Figure 8. The drift of the latitude of maximum (left) and the centroid position of the magnetic flux at the near-surface layer in the 
models D(l-4). The dash-dotted line shows results for the model Dl, the red dashed line - for the model D2, the solid black line - for 
the model D3, the bl ack dashed l ine - f or the model D4 and the solid green line shows the exponential law of the sunspot area centroid 
drift, as suggested by Hathaway (2011). 



latitude as suggested by Hathaway (2011): 



A c (t) = 28°exp 



12* 
90 



where t is time measured in years. Model D4 shows a change between fast (nearly steady dynamo wave) 
drift at the beginning of the cycle to slow drift at the decaying phase of the cycle. The overlap between 
subsequent cycles is growing from model Dl to model D4. In all the models the highest latitude of the 
centroid position of the toroidal magnetic flux is below 30°. Models D3 and D4 have nearly equal starting 
latitude of the centroid position. It is about 24°. This means that a model with increased magnetic helicity 
pumping produces nearly the same effect for the shift of the centroid position as a model with a strong 
nonlinear effect of magnetic helicity. 
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3 Discussion and conclusions 



We have shown that the interaction of helical convective motions and differential rotation in the solar 
convection zone produces a turbulent drift of large-scale magnetic field. The principal direction of the drift 
corresponds to the direction of the large-scale vorticity vector. The large-scale vorticity vector roughly 
follows to iso-surfaces of angular velocity. Since the direction of the drift depends on the sign of helicity, 



the pumping effect is governed by the Parker- Yoshimura rule (Parker 1955, Yoshimura 1975) 



The effect is computed within the framework of mean-field magnetohydrodynamics using the minimal 
r-approximation. In the calculations, we have assumed that the turbulent kinetic and current helicities 



are given. The calculations were done for arbitrary Coriolis number. In agreement with Mitra et al. (2009) 



and [Rogachevskii et al. (2011), the analytical calculations show that the leading effect of pumping is de- 



scribed by a large-scale magnetic drift in the direction of the large-scale vorticity vector and by anisotropic 
pumping which produces a drift of toroidal and poloidal components of the field in opposite directions. 
The component of the drift that is induced by global rotation and helicity (second line in Eq. Q) is rather 
small compared to the main effect. The latter conclusion should be checked separately for a different model 
of background turbulence, taking into account the generation of kinetic helicity due to global rotation and 
stratification in a turbulent medium. 

We have estimated the pumping effect for the solar convection zone and compared it with other turbu- 
lent pumping effects including diamagnetic pumping and turbulent pumping that results from magnetic 



fluctuations in stratified turbulence (Kichatinov 1991, Pipin 2008). The latter is sometimes referred to 
as "density-gradient pumping effect" (Krivodubskij 2004). The diamagnetic pumping is upward in the 



upper part of the convection zone and downward near the bottom. The velocity field of density-gradient 
pumping is more complicated (see Figure 4). However, its major effect is concentrated near the surface. 
Both diamagnetic pumping and density-gradient pumping effects are quenched inversely proportional to 
the Coriolis number ( Kichat inov| 1 1 99 f , |Pipin|[2008 ). The helicity-vorticity pumping effect modifies the 
direction of large-scale magnetic drift at the bottom of the convection zone. This effect was illustrated 
by a dynamo model that shows a dominant poleward branch of the dynamo wave at the bottom of the 
convection zone. 

It is found that the magnetic helicity contribution of the pumping effect can be important for explain- 
ing the fine structure of the sunspot butterfly diagram. In particular, the magnetic helicity contribution 
results in a slow-down of equatorial propagation of the dynamo wave. The slow-down starts just before 



the maximum of the cycle. Observations indicate a similar behavior in sunspot activity (Ternullo 2007 
Hathaway||2011 ). A behavior like this can be seen in flux-transport models as well (Rempel|2006|). For the 



time being it is unclear what are the differences between different dynamo models and how well do they 
reproduce the observations. A more detailed analysis is needed. 
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Appendix A 

To compute S it is convenient to write equations and ^ in Fourier space: 



d 2 \ - i - 
— + rjz \bj = -ziuj (z) Bi + 



(26) 



+ izi [ Ikjz - q)VM) ~ h (z - q)Vi(q) 



dq + Q5j. 



d_ 
dt 



+ vz 2 )ui = fi + & - 2 (flz) (z x u) 



(27) 



- iTTi f (z)zi / ui(z - q)Vf(q) +u f [z- q)Vi(q 



+ -k (z) (z-B), 



dq 
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where the turbulent pressure was excluded from ([2]) by convolution with tensor Kij(z) — Sij — Z{Z^ Sij 
is the Kronecker symbol and z is a unit wave vector. The equations for the second-order moments that 



make contributions to the mean-electro- motive force(MEMF) can be found directly from (26, 27). As the 
preliminary step we write the equations for the second-order products of the fluctuating fields, and make 
the ensemble averaging of them, 



d_ 
dt 



Ui 



{z)bj (z f )j = iz[Bi (ui (z)uj (V)) - 2 (f2z) eu n zi (u n (z)bj(z') 

+ iz[J [(ui (z) k (z' - q)) Vj(q) - (u, (z) b 3 (z' - q)) V^q) 



(28) 



dq 



n if (z) Zl / U (z - q)bj (z'))Vf(q) + (u f (z - q)bj {z'))Vi(q) dq 



dt 



d_ 
dt 



+ -B lZl (h [z)bj (z')) + Th%(z, z') - (rjz' 2 + vz 2 ) (z) bj (z')) , 

(Hi (z)uj [z')) = -2(i?z) eunzi (u n (z)uj(z')) - 2 (i7z') £ji n z[ (ui(z)u n (z')) 

-iir if (z)zi J (iii {z - q)uj (z'))V f (q) + (u f (z - q)uj (z'))Vi(q) dq 

- i-K jf (z')z[ j (ui{z)ui{z - q))V f {q) + I^Ui(z)u f {z - q))Vi{q) dq 
+ 77$ (z, z') - v (z' 2 + z 2 ) {in (z) uj (z')) , 
bi (z) bj {z')) = Th%{z, z') - {r,z' 2 + r,z 2 ) (bt (z) bj (z') 

+ i z[ J [(bi (z) h (z' - q)) Vj(q) - (bi (z) bj (z' - q)) V^q)] dq, 



(29) 



(30) 



where, the terms Th^p v,h ^ involve the third-order moments of fluctuating fields and second-order moments 

of them with the forcing term. Next, we apply the r-approximation, substituting the Th[* ,v,h ^ -terms by 
the corresponding r relaxation terms of the second-order contributions, 



Th^ ) ^-(m i (z)b j {z'))/r c , 



Th 



(v) _^ (m (z) rhj (zQ) - (rhj (z) rhj (zQ) (0) 



b i (z)b j (z'))-(b i (z)b j (z' 



(0) 



(31) 
(32) 

(33) 



where the superscript . . denotes the moments of the background turbulence. Approximating these 
complicated contributions by the simple relaxation terms has to be considered as a questionable assump- 



tion. It involves additional assumptions (see Radler and Rheinhardt 2007), e.g., it is assumed that the 
second-order correlations in Eq. (8) do not vary significantly on the time scale of r c . This assumption is 
consistent with scale separation between the mean and fluctuating quantities in the mean-field magneto 
hydrodynamics. The reader can find a comprehensive discussion of the r -approximation in the above cited 
papers. Furthermore, we restrict ourselves to the high Reynolds numbers limit and discard the microscopic 
diffusion terms. The contributions of the mean magnetic field in the turbulent stresses will be neglected 



because they give the nonlinear terms in the cross helicity tensor. Also, r c is independent on k (cf, |Radler 
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et al.|[2003| |Rogachevskii and Kleeorin"||2003[ |Brandenburg and Subramanian||2005|) and it is independent 



on the mean fields as well. This should be taken into account in considering the nonlinear effects due to 
rotation. Taking all the above assumptions into account, we get the system of equations for the moments 
for the stationary case: 



Ui 0) bj 0') 



= -2(f2z)eu n zi (u n (z)bj(z')\ +iz / l Biu i (z)uj (z')) 



(34) 



+ iz[ 



l 7T, 



: (z) b t (z f - q)) V 3 (q) - {in (z) b 3 (z' - q)) Vi(q)\ dq 
f (z)z t J [(«, (z - q)b 3 (z'))P f (q) + (% (z - q)b 3 (*')) ^t(«Z)] dq 



+ -B l z l (b i (z)b j (z% 



{rht {z)rhj (z')) 



+ 



-2 (£}z)su n zi (u n (z)uj(z f )) - 2 (f2z!) e 3 i n z\ (ui(z)u n (z')) 
'■Ki f {z)zi J (ui [z - q)uj (z'))V f {q) + (u f (z - q)uj (z'))Vi(q 
^iA z ') z 'l J («< ( z ) u l ( z ~ l)) ^f(Q) + ( u i ( z ) Uf (z - q)) V t (q) 

(m i (z)m j (z'))^ 



(35) 



dq 
dq 



bi{z)b 3 {z')) (b i (z)b 3 (z') 



(0) 



(36) 



+ iz[ 



k (z) h (z' - q)) Vj(q) - (k (z) b 3 [z' - q)) V t (q) dq 



To proceed further, we have to introduce some conventions and notations that are widely used in the 
literature. The double Fourier transformation of an ensemble average of two fluctuating quantities, say / 
and g, taken at equal times and at the different positions x, x\ is given by 



</(*MO> 



/ (z) g (z') ) e i ( - z - x+z '- x '^d 3 zd 3 z'. 



(37) 



In the spirit of the general formalism of the two-scale approximation (Roberts and Soward 1975) we 



introduce "fast" and "slow" variables. They are defined by the relative r — x — x 1 and the mean R 
\ (x + x') coordinates, respectively. Then, Eq. (37) can be written in the form 



(f( X )g( X >)) = J J (V ffc + W-fc + 



(38) 



where we have introduced the wave vectors k — i [z — z') and K — z + z' . Then, following to Subramanian 



and Brandenburg ~(^U(J4), we detine the correlation function ot f and g obtained from (38) by integration 
with respect to X, 



* (/, §, k, r) = f(f(k + \k) g ( -k + \k) ) j 



K. 



(39) 
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For further convenience we define the second order correlations of momentum density, magnetic fluctuations 
and the cross-correlations of momentum and magnetic fluctuations via 



Vij (fc, R) = <S>{u h Uj, fc, R), (u 2 ) (R) = J vu (fc, R) d 3 fc, 

ha (fc, R) = fy, k , R ), ( b2 ) ( R ) = J ki (fc, R) d 3 fc, 

(fe, i?) = bj, fc, H), Si (R) = e ijk J x jk (fc, i?) d 3 fc. 



(40) 
(41) 
(42) 



We now return to equations ([34]), (35) and (36). As the first step, we solve these equations about ft (non- 



linear effects of the Coriolis force) and make the Taylor expansion with respect to the "slow" variables 



and take the Fourier transformation, (39), about them. The details of this procedure can be found in 
( |Subramanian and Brandenbur g 2004). In result we get the following equations for the second moments 



^ = -CD?} (B.k) (v fj - ) + Vf}v u k n - D?}v u k l3 + (43) 



,(o) _ Sjf + jjgkpejfp + -ifykikf 
1 + ^ 



D\y = " r " V /,2 " i ^ - 2 (fl.fe) r c 



1 ijnm ~ °in°jm > ^ \£mp°jm > £jmp°m) 

^ij^nm &nmkikj + S{ m k n kj + &njkik m 2(5 n [j(5j] r 



M 



M = 1 + 4^ 



= ^ + T chilVj t l + T c hljV it l + Tcfc/V/, 



(45) 



These equations were solved with respect to the shear tensor, V^j = V^Vi, by means of perturbation 
procedure. One remains to define the spectra of the background turbulence. We will adopt the isotropic 



form of the spectra (Roberts and Soward 1975). Additionally, the background magnetic fluctuations are 



helical while there is no prescribed kinetic helicity in the background turbulence: 



*3 



m = 



TTij (k) 
TTij (k) 



U {k, R) 

ijpKp 8vrfc4 

___ M(k,R) 



E (fc, R) 
8nk 2 

B (k, R 



2/ £<j iryk 



(46) 
(47) 



where, the spectral functions E(k,R),B(k,R),J\f(k,R) define, respectively, the intensity of the velocity 
fluctuations, the intensity of the magnetic fluctuations and amount of current helicity in the background 
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turbulence. They are denned via 

where /ij^ = (t^^-V X t^ )) and = — — — ~- I n nna l results we use the relation between 
intensities of magnetic and kinetic fluctuations which is defined via B(k,R) = EfipE (k, R). The state 
with e — 1 means equipartition between energies of magnetic an kinetic fluctuations in the background 
turbulence. 



(48) 



Appendix B 



& = wn + 391 > ^^F 1 ~ 1456 ^ 2 + V - 3 ^* 2 ~ 355 ) 

/W = _ 1 ( 5 ((10672I2* 2 + 3872) Q* 2 + 337) arctan ( 2 ^*) 
J3 36540J?* 4 V V ' ' 2J2* 

- 320 ((515J?* 2 - 128) Q* 2 - 895) &Ici& ^ Q } + 3 ((2304J2* 2 + 3380) Q* 2 - 48295) 
/W = (3 ((llfT 2 + 8) J2* 2 - 7)) arCt ^ ( ^ ) - (31^ - 21)) 

& ] = l — A ( (74721T 4 + 501612* 2 - 1685) arctan ( 2jy ) 

Jb 55296Z2* 4 V 2/2* 

- 32(79^ 4 + 1410I2* 2 + 4475) arctai ^*) 

(9872J?* 6 + 183632J?* 4 + 238183/2* 2 + 48295) \ 
+ 3 (J?* 2 + 1) (4J2* 2 + 1) J 

f(7 )_ 1 L^,™^ arctan(^) ((l01^ 2 + 20) Q* 2 - 105) \ 

/e -3272^^^ +3Uj ^ " 35j ^ 3(^*2 + 1) ) 

ir 7) = "ccoi^TTI f (20528JT 4 + 16536J2* 2 - 1685) ^^l 2 ^) 
J7 55296Z2* 4 V ; 2Q* 

- 32 (577Z2* 4 + 3660I2* 2 + 4475) &rcta *^ 



(80528J?* 6 + 282512/2* 4 + 258343J?* 2 + 48295) 
+ 3 " (f?* 2 + l) (4/2*2 + 1) " 



fi 



•(7) _ 



(^ 4 + 7) A „ 2 arctan(^) _ (l3/2* 2 + 15) \ 
32i7* 4 \^ j tt* 3(r2* 2 + 1) J 



